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We present a novel approach to the detection of weak magnetic fields that takes advantage of 
recently developed techniques for the coherent control of solid-state electron spin quantum 
bits. Specifically, we investigate a magnetic sensor based on Nitrogen- Vacancy centers in 
room-temperature diamond. We discuss two important applications of this technique: a 
nanoscale magnetometer that could potentially detect precession of single nuclear spins and 
an optical magnetic field imager combining spatial resolution ranging from micrometers to 
millimeters with a sensitivity approaching few femtotesla/Hz^/^. 

The detection of weak magnetic fields with high spatial resolution is an important prob- 
lem in diverse areas ranging from fundamental physics and material science to data storage and 



biomedical science. Over the past few decades, a wide variety of magnetic sensors have been 
developed using approaches including superconducting quantum interference devices (SQUIDs)\ 
the Hall effect in semiconductors^, atomic vapor and BEC-based magnetometry^"^''-^, and mag- 
netic resonance force microscopy ^"^^ In this article we present a novel approach to high spatial 
resolution magnetic field detection, using systems currently explored as quantum bits: isolated 
electronic spins in a solid. We focus on spins associated with Nitrogen- Vacancy (NV) color cen- 
ters in diamond^^ (Fig. 1 a-b), since they can be individually addressed, optically polarized and 
detected, and exhibit excellent coherence properties even at room temperature^^"^^. Recently, co- 
herent control of NV electronic spin qubits has been used to sense and manipulate nearby individ- 
ual electronic ^ ^' " and nuclear spins^^ in a diamond lattice. Here we describe how such a system 
can also be used for the precision sensing and imaging of external magnetic fields. 

We discuss two types of potential implementations of such sensors. First, a single sensing 
spin confined in a nanoscale region can be brought in direct proximity to a magnetic field source, 
such as an electron or nuclear spin. For example, a diamond nanocrystal (10-50 nm in size) con- 
taining a single NV center can be attached to a tip of a scanning probe (Fig. Ic)^^. Second, a bulk 
diamond sample with a high density of NV centers can be used to sense fields created by remote 
objects with ultra-high sensitivity and sub-//m spatial resolution (Fig. Id). 
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Magnetometry with single electronic spin qubits 



The operating principles of our approach are closely related to those of magnetometers based on 
spin precession in atomic vapors. In particular, detecting the relative energy shift induced by a 
magnetic field b between two Zeeman sublevels allows for a precise determination of an applied 
DC or AC magnetic field. Ultimately, sensitivity is determined by the spin coherence time and 
by the spin projection-noise. Although solid-state electronic spins have shorter coherence times 
than gaseous atoms, quantum control techniques can decouple them from the local environment 
and from each other, as we show below, leading to a substantial improvement in their sensitivity 
to external, time-varying magnetic fields, while retaining the desirable features of a robust solid 
sensor. 

The canonical approach to detecting a Zeeman shift uses a Ramsey-type sequence as il- 
lustrated in Fig. |2^. A 7r/2-pulse creates a superposition of two Zeeman levels, which acquire a 
relative phase (f) = 6lj r oc ^jf-br from the external field b during the free evolution interval r (here 
/iB is the Bohr magneton and (7 ^ 2 for NV centers). Another 7r/2-pulse transforms the relative 
phase into a population difference, which is measured optically and from which the Zeeman shift 
is inferred. For small (p, the magnetometer signal S (proportional to the induced population differ- 
ence) depends linearly on the magnetic field: S ~ ^^br. During the total averaging interval T, 
T/r measurements can be made, yielding a shot-noise-limited sensitivity rj given by the minimum 
detectable field, bmin = v/VT = 

Increasing the interrogation time r improves the sensitivity until random (environmental) 
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perturbations lead to decay of the free-precession signal. In the case of solid-state spin systems, 
the coherence is limited by interactions with nearby lattice nuclei and paramagnetic impurities, 
resulting in an ensemble dephasing time T2*. Furthermore, there will be a finite number of fluores- 
cence photons collected and detected, leading to additional photon shot noise, and a finite contrast 
to the Ramsey fringes. We describe these effects by a single parameter C < 1, which approaches 
unity for ideal, single-shot readout (see Methods). The optimum sensitivity of a magnetometer 
based on a single electronic spin, achieved for r ~ T2 , is given by 



For current experiments^^, with detection efficiency ~ 10"^, C ^ 0.05 and ~ 1 /xs. This 
yields an optimal sensitivity ~ 1 //T/Hz^/^. Improving the collection efficiency to 77 ~ 5% yields 
C !^ 0.3 and leads to a sensitivity ~ 120 nT/HzV2. 

Coherent control techniques can improve the sensitivity for AC fields. Due to the long cor- 
relation times characteristic of dipolar interactions between spins in systems such as diamond — 
the principal source of dephasing — spin echo techniques can dramatically extend the coherence 
time. Specifically, by adding an additional microwave tt pulse to the Ramsey sequence at time 
t/2, the Hahn echo sequence (Fig. 2a) removes the effect of environmental perturbations whose 
correlation time is long compared to r. Thus a signal field b{t) oscillating in-phase with the 
pulse sequence produces an overall additive phase shift, leading to a total phase accumulation, 
Scf) = sm[J^/^ b{t)dt - jl^^h{t)dt]. For a signal field h{t) = bsm{ut + (po), this yields d(f) = 
SJfbrfiur, ifio), with f{x, ipo) = «i'^'(^/^)^°«(^/2+yo) ^ essence, the spin echo allows us to extend 
the interrogation time r from the limit set by T2 up to a value T2 which is close to the intrinsic 
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spin coherence time, at the cost of a reduced bandwidth and insensitivity to frequencies < I/I2. 
For maximal response to CW signals with known frequency and phase (assuming small b), we find 
T = 27: /u and </7o = to be optimal. For signals whose time dependence is a-priori unknown, it 
is useful to measure the signal variance, which provides equivalent performance (see Methods). In 
either case, the sensitivity is improved by ~ ^/TfJT^^. 



The optimum sensitivity is achieved only for fields oscillating near 1/ ~ 1/72. However, 
these results can be easily extended to higher frequency signals. In particular, for signal field 
oscillation periods shorter than the dephasing time, the interrogation time need not be restricted 
to the duration of one period, but can be multiples of it. Then, composite pulse sequences such 
as CPMG^^ may perform better at the expense of a reduced bandwidth. Furthermore, in ultra- 
pure samples where nuclear spins' evolution leads to decay of the echo signal, the long correlation 
time of the nuclei leads to non-exponential decay of the echo signal^^'^^. In this case, the CPMG 
sequence can increase the interrogation time, further reducing the minimum detectable field (see 
Fig. 2 and Methods). Finally, another way to improve the magnetometer sensitivity is to use many 
sensing spins, where we can take advantage of the relatively high achievable density of spins in the 
solid-state (~ 10^^ cm~^) compared to atomic magnetometers(~ 10^^ cm"^).^"^ 



(2) 
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Implementation with NV centers 

We now discuss specific details of magnetometry using Nitrogen- Vacancy centers in diamond, 
developing an appropriate operating regime and determining the optimal sensitivities possible for 
current experimental technology. The fine structure of the electronic ground state of a NV center, 
shown in Fig. [T]-a, is a spin triplet. The crystal field splits the m., = ±1 Zeeman sublevels from 
the uis = sublevel by A = 27r x 2.87 GHz, allowing the use of electron-spin resonance (ESR) 
techniques even at vanishing external magnetic field. Furthermore, under application of green light, 
NV centers exhibit a transient, spin-dependent fluorescence, which allows for optical detection of 
the spin. After the transient signal decays, the system optically pumps into the = state, which 
prepares the system for the next measurement (see Methods). 

As a specific example, we focus on magnetometry in low external static magnetic fields (< 
10 mT). In this case, A is the largest energy scale and sets the spin quantization axis parallel 
to the nitrogen-to-vacancy direction. The secular Hamiltonian, including a small external field 

B{t) = {B,, By, B,), is 

n = hAS',+gfiBB,S,, (3) 

where B^ is the component of the magnetic field along the NV center's axis and takes the values 
rus = 0, ±1. Terms proportional to the perpendicular field are suppressed to order ~ B^ y/ A and 
do not depend on the field B^ being measured, and therefore may be neglected. 

At low magnetic fields, the = ±1 manifold can be used to implement a vector magne- 
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tometer, sensitive only to components of the magnetic field along the center's axis. We focus on 
the nig = ±1 manifold as it has twice the energy splitting of the 0-1 manifold and is less affected 
by nuclear spin-induced decoherence at low fields, since inter-nuclear interactions are suppressed 
by the large hyperfine field^^. 

Coherent control of the NV center's spin states is obtained via an ESR magnetic field os- 
cillating at angular frequency A. ESR pulses linearly polarized along the x-axis rotate the NV 
spin between the two dimensional subspace of |0) and |+) = (|1) + | — l))/v^. To manipulate 
I ±1 ) superpositions, additional control can be provided by a background oscillating reference field 
[Srcf sin(27rt/r)] along the 2;-axis. Specifically, i?rcf = {fi/ QI^b)'^'^ /^t yields an optimal phase 
offset to achieve a magnetometer signal linear in the field strength (Fig. 2). 

The sensitivity as a function of the signal frequency for both AC and DC detection is plotted 
in Fig.[2l For diamond where natural abundance (1.1%) Carbon-13 nuclei are the principal cause 
of dephasing, Tg* ~ 1/iS and T2 ~ 300/is'^. Again using current experimental parameters, with 
C ~ 0.05, and measurement and preparation time tm < 2 yus, we can optimize the sensitivity as a 
function of r. Including corrections from decoherence with expected signal decay (see Methods) 
cx exp[— (r/T2)^] we find: tiac = ^^e^^/^^)'^^-^ _^ tm/Cr. We obtain optimal sensitivity of 
r]Ac ~ 18 nT Hz^^/^ for a single NV center using current experimental collection efficiencies. 
Improved collection efficiencies (C = 0.3) would yield -qAc = 3 nT/Hz"^/^. Note that spin Ti 
relaxation occurs on timescales much longer than milliseconds and may be safely neglected^^. 
Finally, the observed dephasing times are independent of temperature from 4 K to 300 K, due in 
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part to the vanishing polarization of the nuclear bath at small magnetic fields. 

When more than one nitrogen vacancy center exists in the sample, they can belong to four 
different crystallographic classes, each corresponding to the centers' alignments along different 
(111) axes. To operate as a vector magnetometer along a controlled direction, a transverse (DC) 
magnetic field B± > 0.3 mT (see Methods) detunes the other three classes' levels such that the 
ESR field used for quantum control excites only spins with the desired crystallographic orientation, 
perpendicular to the external field. Thus 1 in 4 spins contribute to the magnetometer signal. 

Magnetometry in the high density limit 

A principle advantage of our approach over other spin precession magnetometers is the high 
achievable density n of sensing spins. This improves the sensitivity to fields that are homogeneous 
over the magnetometer volume, since the projection-noise per unit volume decreases as \l\fn. 
NV centers can be created in controlled densities by implanting high-purity diamond with nitrogen 
ions and subsequently annealing the sample to recombine the nitrogen with vacancies^^. Assuming 
an initial nitrogen concentration ~ 10^^ cm~^ with a conversion / ~ 0.1 to NV centers we 
expect it will be feasible to create diamond crystals with an NV center density exceeding ~ 10^^ 
cm~^, with an average distance between centers of less than 10 nm. Even at these densities, effects 
such as superradiance do not play a role due to the large spectral width of the NV fluorescence. 

At high spin densities, NV-paramagnetic impurities and NV-NV interactions may limit the 
sensitivity of the magnetometer. In particular, substitutional Nitrogen impurities with one bound 
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electron (PI centers) become a sizable source of dephasing in high density samples^^'^^. The 
dipole-dipole interaction between these centers has a characteristic time scale = — , where 
a is on the order of the dipole coupling between electron spins, ^ '■^^^^ ^ 3.3 x 10"^^ s^^cm^, 
and n^pr is the density of paramagnetic impurities. Qualitatively, this time scale corresponds to 
the rotation time of a single paramagnetic spin in the presence of the random field from the other 
paramagnetic centers. The time scale for interaction between this impurity bath and a given NV 
center will be of the same order of magnitude. This suggests an exponential decay of spin echo 
coherence on a timescale (see Methods), in contrast to single NV center-based sensing where 
nearby nuclear spins limit the coherence time. 

To evaluate the effects of paramagnetic impurities, we assume a density n of NV centers and 
Uepr = n{\ — /)// of paramagnetic impurities, where / is the conversion factor described above. 
The relevant figure of merit is the sensitivity per root volume rj\(j = tjacW. We find 

v\c = rr- ^ (4) 

gjj^B Cy'riT 

where we have taken into account that the sensing centers account for only one fourth of the 
NV centers in the sample. Here we include both dephasing due to a bath of dipolar-coupled 
nuclear spins and the paramagnetic spin bath just discussed. In the high NV density- and low 
/-regime, T2^carbon > Tc > T^, i.e.. Carbon- 13 is no longer the limit to echo lifetimes, but 
still limits inhomogeneous broadening. Then the optimum magnetometer sensitivity becomes: 



ri\^* = -I^^^£^<tJl, For / = 0.1 and T2,carbon = 300 yus, the optimum sensitivity is NV 
density independent over the range n ~ 10^^ — 10^^ cm^^, as is seen in Fig. [3^, and reaches a 
maximum sensitivity value r]v ~ 250 aT Hz^^/^ cm^'^/^ for C = 0.3. However, the optimum echo 



time depends upon the NV density, r = //[(I — f)2na], with higher density samples requiring 
higher detection frequencies. Finally, for n ^ 10^^ cm^^, corrections due to finite preparation, 
control, and measurement times can become important, and lead to the limitations in sensitivity at 
high NV density seen in Fig. [3^. 

To push the sensitivity limits beyond the cutoff imposed by paramagnetic impurities, we can 
exploit more advanced forms of dynamical decoupling^" than spin echo. With appropriate exter- 
nal time-dependent controls, the system can be made to evolve under an effective, time-averaged 
Hamiltonian that is a suitable symmetrization of the undesired interactions. For example, driving 
the PI centers via spin resonance at a rate much faster than the intrinsic decorrelation time, Tc, 
acts as a rapid spin-echo for the NV centers without impacting the NV center's magnetic field 
sensing capabilities. Furthermore, improving implantation and conversion techniques (by optimiz- 
ing implant energies^' or by using cold implantation^') could increase the ratio of NV centers to 
paramagnetic impurities. When the conversion efficiency exceeds 50%, interactions between NV 
centers become the primary source of noise, with a dephasing^^ oc (anr)^. The coupling between 
the sensing NV centers is a Sj^zSk,z interaction which is not removed with spin echo. However, 
by using collective rotations driven by appropriate ESR pulses, the interaction can be successively 
rotated through the x, y and z axes for an equal time duration^^; so that on average the spins will 
experience an isotropic Hamiltonian, which commutes with the signal perturbation and thus allows 
the spin evolution necessary for magnetometry^'^. Pulse sequences such as MREV^^'^^ can induce 
the desired evolution, and will be necessary in the high NV-center density limit. 
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Single spin detection with a nano-NV-magnetometer 

NV magnetometers can be applied to an outstanding challenge in magnetic sensing: the detection 
and real space imaging of small ensembles of electronic and nuclear spins, with the long-term goal 
of resolving individual nuclear spins in a molecule. Since the magnetic field from a single dipole 
decreases with distance as ~ a magnetometer that can be brought into close proximity of the 
field source offers a clear advantage. A diamond nanocrystal or a single NV center near the surface 
of a bulk crystal would allow for a spatial resolution limited only by the distance between the NV 
center and the object of study, not by the wavelength of the fluorescence signal. For example, 
consider as a prototype system consisting of a crystal with a single NV center at a distance tq ~ 10 
nm from the surface of the crystal. At this distance, the dipolar field from a single proton is i?H — 3 
nT, which is well within the projected limits for a single NV center. 

To examine a practical method to measure the magnetic field from a single spin, we consider 
a material with a varying nuclear spin density Ug that is brought in close proximity (a distance 
~ To) to the NV center. At realistic temperatures, the thermal nuclear spin polarization of the 
material will be small. However, because only a few spins are involved, the distribution of spin 
configurations leads a large variance in the spin polarization^^, providing a substantial, albeit ran- 
domly oriented, magnetic field to be detected by the NV magnetometer. We find (see Methods) that 
the field magnitude measured by our sensor will be characterized by a variance i?rms ~ B^iy/N, 
where N ~ STr^sro is the effective number of spins contributing to the signal. This indicates that 
our prototype system has an effective spatial resolution determined only by the distance of the NV 
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center from the surface of the sample material, assuming one can position the sensor relative to the 
sample with stability much better than tq. 

At nuclear spin densities < 10^*^ cm~^, there is on average one or fewer nuclear spins in an 
effective sensing volume with ro ~ 10 nm. Hence, in this case single spins could be measured. 
However, most organic molecules have substantially higher proton densities (~ 10^^ — 10^^ cm^^). 
To measure only one proton at a time would require a further improvement in spatial resolution. In 
this case a magnetic field gradient can be used, which allows one to convert high spectral resolu- 
tion to high spatial resolution. Using similar techniques to those used in magnetic force resonance 
microscopy^, a magnet near the surface of a substrate can produce gradient fields of > 10^ T/m 
(Fig. lib). The narrow bandwidth of our detector, ~ 4 kHz (~ 1/^2)5 allows it to spectrally dis- 
tinguish two protons separated by a magnetic field difference of 0.1 mT, corresponding to physical 
separation of 0.1 nm. This implies that individual proton detection may be possible even in organic 
and biological molecules. The narrow bandwidth associated in particular with the CPMG approach 
(see Methods) allows one to distinguish different isotopes, due to their unique gyromagnetic ra- 
tios. More generally, our approach enables the detection of nanoscale variations in the chemical 
and physical environment. 

We note that the present approach can surpass the sensitivity of SQUID\ Hall-bar^, and re- 
cently proposed optically-pumped semiconductor-based^^ nano-magnetometers by more than an 
order of magnitude, with 10-1000 times better spatial resolution. The ultimate limits to miniatur- 
ization of NV center nano-magnetometers, which are likely due to surface effects, are not yet well 
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understood. 



Imaging of macroscopic magnetic fields 

In contrast to the nano-magnetometer approach outlined above, a macroscopic crystal of diamond 
containing many NV centers may be used as a high sensitivity imaging magnetometer with large 
field-of-view and optical wavelength-limited spatial resolution. As an example system we consider 
a crystal of diamond with a high density of NV centers. The signal from NV centers in a diffraction 
limited setting, where a CCD might be used to image the crystal, is divided into separate "pixels", 
with each pixel corresponding to a ~ (1/xm)^ volume element of the crystal. For NV center 
densities of ~ 10^^ - 10^^ cm-=^ and C = 0.3, each pixel would have ~ 100 pT Hz-V2 AC 
sensitivity. This spatial resolution is comparable to micro-SQUID magnetometers but with four 
orders of magnitude higher magnetic field sensitivity In such a scenario, diamond crystals 
could range from tens of microns to millimeters in size; and be physically integrated with fiber- 
based optics for a robust and practical magnetic field imager. 

Larger detector volumes further improve the sensitivity for whole-sample measurements. For 
example, a (3 mm)^ x 1 mm thick crystal can achieve an overall sensitivity of 3 fT Hz~^/^ with mm 
resolution. Reducing the ratio of paramagnetic impurities to NV centers could potentially lead to 
the detection of attotesla fields, opening the prospect of improved tests of fundamental symmetries 
and physical laws. 
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Conclusions 



The extremely high magnetic field sensitivity in a small volume offered by solid state spin-qubits 
such as NV centers in diamond can find a wide range of applications, from fundamental physics 
tests or quantum computing applications to detection of NMR signals, surface physics and mate- 
rial science, and medical imaging and biomagnetism. Recently, a proof-of-principle experimental 
demonstrations of such a sensor have been performed by members of our collaboration''^^ and other 
groups'*^ . Further extensions could include the use of non-classical spin states, such as squeezed 
states induced by the spin-spin coupling. The sensitivity could also be improved by using synthe- 
sized, isotopically purified diamond containing a lower fraction of Carbon- 13, the main cause of 
dephasing at moderate NV densities, and by developing more efficient NV center creation tech- 
niques that do not result in high densities of paramagnetic impurities. On a more general level, 
these ideas could apply to a variety of paramagnetic systems or other types of solid-state qubits 
that are sensitive to different perturbations. 

Methods 

ESR control techniques The NV center's spin triplet has a V-type level configuration. An external 
microwave field tuned to the A = 2.87 GHz resonance with its magnetic field linearly polarized 
along the x-axis drives transitions between |0) and the superposition |+) — + |— l))/v^, 
while the state |— ) = (|1) — \—l))/\/2 is dark — it is decoupled from the field due to destructive 
quantum interference. Application of a magnetic field aligned with the NV-center 2; -axis perturbs 
the interference, and allows for complete quantum control of the spin triplet. In an echo sequence 
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appropriate for magnetometry using the |+) and | — ) states, the traditional 7r/2 — tt — 7r/2 structure 
is replaced by tt — 27r — tt: the first pulse creates |+), the second induces a relative tt -phase shift 
between |+) and |— ), and the third converts |+) to |0) while leaving |— ) population trapped in 
the rus — ±1 manifold. We remark that for external fields in excess of a few mT it may be more 
convenient to use the 0- 1 manifold, as two different resonance frequencies would be necessary for 
using the ±1 manifold in this regime. 

AC-field measurement scheme and bandwidth AC-field detection requires synchronization of 
the pulse sequence with the external magnetic field oscillations. When this is not practical or if 
the field phase ^po varies randomly in time, successive measurements will give random readings 
distributed over the range of the function /{ur, (po) (given in the main text) resulting in a zero 
average signal. In this situation, information about the field intensity is contained in the measured 
signal variance, provided the random phase correlation time satisfies: r <^ < T. (If > T, 
the total averaging time, the scheme presented in the main text could be used). For r = 27r/u 
and a uniformly distributed (po, (/(27r, <^o)^) = 2/7r^ and the standard deviation of the measured 
signal is: ^^^^ br, while the noise has a contribution from the uncertainty in the variance equal 



to 2^/^/7r. The sensitivity is thus only worsened by a factor y 2(1 + v^/tt^) 1.5 compared to 
detection of a signal with a known phase. 

To increase the sensitivity at higher frequencies, it is possible to increase the interrogation 
time (see main text) by using a series of 27r-pulse cycles (CPMG pulse sequence). A single cy- 
cle corresponds to the pulse sequence t/4 — tt — t/2 — tt — r/4. While this method increases 
the sensitivity, the measurement bandwidth decreases with increasing cycle number Uc- The AC 
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magnetometer response to a general signal b{t) can be calculated from a frequency space anal- 
ysis: ^ Too kuj) e-*dt - e-*dt) da; = ^ JZ, b{uj)Woiu;, T)e^-^lHu, where is 
the Fourier transform of the signal field and Wq a windowing function. With a similar calculation 
we obtain the windowing function for an ^c-cycle pulse sequence: 

TLU/Z 

This function has a band-center k,2tx/t and bandwidth (HWHM) ~ 4/ {uct). 

We can evaluate the improvement in coherence times for the CPMG sequence in cases 
where a detailed understanding of the main source of decoherence is available. For the single- 
spin magnetometer, we can approximate the nuclear spin environment by separating contribu- 
tions from distant nuclear spins, undergoing dipolar spin diffusion, and nearby nuclear spins, 
whose evolution is frozen by the electron spin's dipolar field^^-^^. We can model the distant nu- 
clear spins as an exponentially correlated gaussian fluctuating field B with a correlation function 
{B{t)B{t')) ~ {-^j^f exp(-|t - t'\/Tc) where > is the correlation time of the nuclear 
spins. 

Within this model, the random phase accumulated during an echo sequence (S(f) — J^^^ B (t) dt— 
j^j^ B{t)dt) is characterized by its variance, ((50^) fa /[QTc{T*f] for > T*, r. Applying this 
model to an nc-cycle CPMG sequence gives (50^) ~ (ncr)^/[24n^Tc(T2*)^]. Thus the multiple- 
pulse sequence yields an improvement in the lifetime by {2ncY/^?^ Th^ improvement is condi- 
tional on r, T2 <C Tc and on the total interrogation time UcT being less than the relaxation time of 
the electron spins. Recent experiments have shown that the relaxation time in ultra-pure samples 
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is ^ 20 ms^^, suggesting ric > 40 cycles can result in an (2nc)^/^ > 4 overall improvement in 
sensitivity. Note that in practice this improvement will be limited by imperfections in the con- 
trol pulses. For example, vr -pulse errors of order 1% will limit Uc ~ 25, resulting in the optimal 
sensitivity show in Fig.|2l 

Measurement efficiency The state of the electronic spin is measured by spin- selective fluores- 
cence. When illuminated by green light, NV centers in the = state undergo a cyclic 
transition'^^, with a rate limited by radiative decay (7 ~ 15 MHz). At the same time, centers 
in the = ±1 state are rapidly pumped into a dark singlet state, from which they decay to the 
nis = state after a time tp ^ 0.5 fis. To allow for a good discrimination of the = 0, ±1 
states, the measurement time tm should be smaller than the optical pumping time tp. 

For a given photon collection efficiency r]m, an average of «o — {tml)Vm photons are de- 
tected from each spin in the = state and ai(< ao) photons are detected from each spin in 
the rris = ±1 manifold. We can estimate the combined effects of spin projection-noise and photon 
shot noise for measurements as N^^^"^ /C, recovering the formulae for sensitivity used in the 



main text, with 1/C = + 2 , . This includes the effects of photon shot noise and reduced 



contrast. For current experiments, a contrast (ao — a;i)/(a;o + «i) ~ 0.3 is observed. Efficiencies of 
r]m ~ 0.001 are achieved in current experiments'^' and give C ~ 0.05. Assuming high collection 
efficiency (77^ > 0.05) gives C ~ 0.3. 

Effects of different NV center orientations In order to use an ensemble of NV centers as a vector 
magnetometer, the signal should originate only from one of the four different crystallographic 




17 



axes. Under application of a DC transverse magnetic field B±x, the other (spectator) centers 
(with crystalline axis h) have their |±1) levels split by qubBx^x ■ n. This detunes the spectator 
centers from the microwave field used for preparing and manipulating the ms = ±1 subspace. For 
example, to use NV centers along the (1,1,1) crystallographic axis, the ideal choice of x is to align 
it with the (1, 1, 2) axis. We require the microwave Rabi frequency Vt > Stt/Tj for pulse errors 
to be smaller than our assumed measurement errors for the desired (111) axis. This translates 
to a requirement that giJiBB± > 3h^^/3/2 for the other three axes. For T2 = 1/xs, we require 
B± > 0.3 mT. One intriguing development of NV center based magnetometry would be to exploit 
the four crystallographic classes of NV centers to provide a full (3D) vector magnetometer, by 
changing the direction of the biasing transverse field B± in between measurements. 

Errors due to inhomogeneities in the NV center properties (e.g., variations of the -factor due 
to crystal strains) or to spatial inhomogeneities of the magnetic field can typically be neglected. 
Even for an average microtesla signal field, a distribution of gi-factors of field inhomogeneity of 
order induced dephasing is 4% leads to a broadening of the signal that is smaller than the effects 
ofTs. 

Coupling to paramagnetic impurities The coupling of an NV electronic spin to other NV centers 

— * — * 

(Sk) and paramagnetic (epr) impurities such as nitrogen (Ik, gi ~ g) is given by the magnetic 
dipolar interaction. To first order in 1/ A, the secular dipolar Hamiltonian is given by: Hzz + ^epr> 
with Hzz = Y,jk i^zjDjk ■ ZkSz,k and Hepr = Y.jk ^zjDjk ■ Ik- The dipole interaction vector is 
^jk — '^°lJn^ {^{■^3k-z)rjk-z\ ^ |.]jg - gg|. |.]jg ]sj_Y crystal axis of the sensing spin centers. 



18 



We model the secular component of the dipole coupling between paramagnetic impurities 
and NV centers as Ujk = Dj^ ■ xl^ to the jth NV center, and with a characteristic correlation time 

can now calculate the expected 

spin-echo signal as a function of (-D^), which scales as the square of the density of paramagnetic 
impurities. In this limit, when the correlation time and the interaction energy are at comparable 
scales, spin echoes decay exponentially as exp(—t/Tc). We find ^ A/Vom^; hence for param- 
agnetic impurity densities of 10^^ cm~^, Tc ^ 1 /is. 

At high densities, paramagnetic impurities and spectator NV centers may have sufficiently 
strong interactions to reduce the correlation time of the field- aligned component, I^- Spectator NV 
centers may be optically pumped to their nis — state, reducing dynamical noise reducing the 
effective temperature of the spectator system. However, spin echoes will not remove the effects 
of the paramagnetic impurities with short correlation times, and they may in fact limit the T2 
time and the corresponding bandwidth of the system. Experiments in systems with high nitrogen 
concentrations indicate exponential decay of echoes on a 5 — 10 /xs timescale^*^ '*^ due to this 
coupling; more generally, the decay scales with the density of impurities. While approaches such 
as CPMG and more complex decoupling may help, we anticipate that paramagnetic impurities 
concentrations below 10^^ cm~^ will be necessary to achieve the best predicted sensitivities of this 
paper. 

Number of spins detected by a point-like sensor To estimate the number of spins that a local- 
ized sensor will detect we determine the maximum and root-mean-square magnetic fields from 
a randomly distributed set of dipolar spins. We denote the dipolar field at a position ro from 
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a spin i as Qb{fi — tq, /,), with normalization of Q = = Bii(10 nm)^ for protons and 

b{r, I) = ^(-T— 3r(r • I)/r'^) being the position dependence of the dipolar field (i?H is the magni- 
tude of magnetic field created by a proton at a distance of 10 nm). The maximum detectable field 
occurs for polarized spins pointing perpendicular {z axis) to the surface (Fig. |4^). By symmetry, 
this field is parallel to the polarization, and we find 




pos 



= -2'nGIn 



where we choose coordinates such that the half-plane begins at z = —tq, ( )pos averages over a 
homogenous distribution of spin positions, and we assume a density ng of dipolar spins, allowing 
us to replace the sum Yli with an integral Ug J^^^^^ d^r. At high temperatures, the fluctuations of 
the potential values of the dipolar field reflect the noise statistics from a set of N spins. The 
mean-square of the 2;-component of the magnetic field is then: 




where the average over spin configurations at high temperature uses (/i,^/j>)c/c/ = ^fiu^ij 3 - 
We find in particular that the statistical fluctuations are consistent with i?rms ~ BhVN, where 
N ~ rigrQ. More specifically, the effective number of spins detected can be estimated from the 
relation B,^^ = \B^^^\/^/N. Thus, = {\B^^^\/ B^msf = j(4i)/3 (87rn,rg). For / = 1/2 this 
reduces to A^ = STrn^rg, equivalent in effective detection volume to a half-sphere of radius 2.29ro. 
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Figure 1 : Overview of a diamond-based magnetometer, (a) Level structure of a single NV center. 
The NV center's ground state is a spin triplet with a A = 2.87 GHz crystal field splitting and a 
Zeeman shift 5uj. Under the application of green light (~ 532 nm), the NV center initially exhibits 
spin-dependent photo-luminescence, even at room temperature, allowing for optical detection of 
electronic spin resonance. After continued illumination the NV spin is pumped into the ground 
state TUs = 0. (b) Crystal structure of diamond with a (111) NV center. A static bias field i?_L is 
applied perpendicular to the 111 axis, while small magnetic fields aligned with the 111 axis are 
detected as the signal, (c) A nano-crystal of diamond at the end of a waveguide for photon collec- 
tion, with resolution limited by the size of the crystal, or (d) a macroscopic sample of diamond, 
with resolution limited by optics, allows for high spatial resolution and signal-to-noise. A green 
laser produces spin-dependent photo luminescence, detected by measuring red light imaged onto a 
CCD. 
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Figure 2: a) ESR pulse sequences for magnetometry, where x and y indicate the linear polarization 
of the ESR pulse in the lab frame. Left: Ramsey pulse sequence for DC-field measurement. Mid- 
dle: echo-based pulse sequence for AC magnetometry 7r/2|^ — n\x — n/2\y. Right: CPMG-based 
pulse sequence for improved AC magnetometry 7r/2|a;(-7r|x- — 'n-\x-)""'7r /2\y, where Uc is number 
of repetitions of the paired vr pulses. For small accumulated phases, a signal linear in the field 
can also be obtained with all pulses along the x direction if a reference field i?ref sin(27rt/r) is 
added, b) DC and AC sensitivity to magnetic fields for a single NV center as a function of signal 
frequency, v. Also shown is the expected performance of CPMG composite pulse sequences, with 
the optimum Uc as described in Methods. Parameters used assume Carbon- 13 limited coherence 
with T* = 1 /is,i3 ^ ^s,^5 ^^^^^^^q ^ o.3, Ti = 20 ms,^^ and an error per pulse of 1%. 
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Figure 3: Sensitivity-per-root-volume (1]^^) at high NV center density, for the AC-field echo mea- 
surement scheme. Left: Contour plot of log^g r/^ in T Hz^^/^cm'^/^ as a function of NV center 
density and signal field frequency. Right: sensitivity at the optimal field frequency, as a function 
of NV center density; the black curve is the sensitivity for / = .1 while the blue and red curve 
are for / = .05 (higher paramagnetic impurity density) and .5 (lower paramagnetic impurity den- 
sity), respectively. Parameters used correspond to T2^carbon = 300 /is,'^, tm = I yus, C = 0.3 and 
a = 3.3 X 10-^3 s-^cm^. 
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Figure 4: Illustration of high spatial resolution magnetometry with a diamond nanocrystal. (a) The 
dipolar fields from spins in the sample decay rapidly with distance; only those within a distance 
~ To contribute to the observable signal for a point-like detector (such as a single NV center in a 
nanocrystal, illustrated by the blue dot). The inset shows how -Bmax and B^^^s are related; when 
few spins are involved, the statistical fluctuations become large, (b) In the presence of a magnetic 
field gradient (field lines in gray) only a small region of the detection volume is precessing at the 
frequency band-center of the detector, allowing for even higher spatial resolution. 
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